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ABSTRACT 


Although sea slaters Ligia have a significant role in rocky shore habitats, their taxonomic entities have not been clearly 
understood. In this study, we investigated whether genetic variation inferred from a nuclear genetic marker, namely 
amplified fragment length polymorphism (AFLP), would conform to that of a mitochondrial DNA marker. Using 
both the mitochondrial DNA marker and the AFLP marker amplified by the six selective primer sets, we analyzed 
95 Ligia individuals from eight locations from East Asia. The direct sequencing of mitochondrial 16S rRNA gene 
revealed three distinct genetic lineages, with 9.8-11.7 Kimura 2-parameter genetic distance. However, the results of 
AFLP genotyping analysis with 691 loci did not support those of mitochondrial DNA, and revealed an unexpectedly 
high proportion of shared polymorphisms among lineages. The inconsistency between the two different genetic 
markers may be explained by difference in DNA evolutionary history, for example inheritance patterns, effective 
population size, and mutation rate. The other factor is a possible genomic island of speciation, in that most of the 
genomic parts are shared among lineages, and only a few genomic regions have diverged. 
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INTRODUCTION 


Sea slaters Ligia are widely distributed in rocky intertidal 
shores, and are significant in the study of the evolutionary 
transition of animals that have successfully adapted from 
aquatic to terrestrial environments. Several aspects of the 
life history of sea slaters suggest that their dispersal capa- 
bility is probably extremely limited, which leads to the pre- 
diction that substantial genetic differentiation should occur 
throughout their range. In these matters, many phylogenet- 
ic studies have recently been conducted on various Ligia 
species worldwide (Jung et al., 2008; Hurtado et al., 2010, 
2018; Markow and Pfeiler, 2010; Eberl et al., 2013; Yin et 
al., 2013). 

Twelve species of Ligia has been recorded in East Asia, so 
far (Nunomura, 1979, 1983, 1990, 1999; Lee, 1994; Schmal- 
fuss, 2003; Tsuge, 2008; Via and West, 2008; Nunomura et 
al., 2011; Via, 2012) and most of these species are allopat- 
rically distributed. However, the two species, Ligia exotica 


Roux, 1828 and Ligia cinerascences Budde-Lund, 1885, are 
sympatrically distributed across East Asia, from the southern 
part of China, to Korea, Japan and Russia. Notwithstanding 
their wide distribution and evolutionary importance, there is 
still uncertainty about the identity of species composition, 
which cannot be resolved by traditional approaches, due to 
their great morphological variation, and large number of in- 
dividuals (Jung et al., 2008). 

Until recently, sequence analyses of mitochondrial DNA 
genes have been the most common method for the phyloge- 
ny and species delimitations of Ligia, especially using COI 
gene and 16S rRNA gene (Jung et al., 2008; Hurtado et al., 
2010, 2018; Markow and Pfeiler, 2010; Eberl et al., 2013; 
Santamaria et al., 2013; Yin et al., 2013). Previous studies 
have proposed distinct genetic lineages in East Asia (Jung et 
al., 2008; Yin et al., 2013). 

Usually, studies using only mitochondrial DNA may cause 
biased results because of its maternal inheritance patterns. 
Thus, an examination of both mitochondrial and nuclear 
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Table 1. Sampled Ligia specimens with geographic coordinates and number of individuals of which 16S rRNA gene and AFLP are 


analyzed in this study 


No. Sampling sites Geographic coordinates Lineages 16S rRNA? AFLP? 
South Korea 
1 Eulwang WG_EW 37°26'37"N, 126°22'38"E WG 11 11 
Uljin EG_UJ 36°59'15"N, 129°25'04"E EG 12 12 
3 Jeju EG_JJ 33°13'04"N, 126°30'49"E EG 12 12 
Japan 
4 Niigata EG_JNI 37°56'04"N, 139°03'08"E EG 12 12 
5 Sendai WG_JSN 38°19'54"N, 141°03'08"E WG 12 12 
China 
6 Qingdao EG_QD 36°05'25"N, 120°29'21"E EG 12 12 
7 Zhujiazien EX_ZH 29°52'12"N, 122°23'55"E EX 12 12 
Russia 
8 Vladivostok WG_VB 43°04'27"N, 131°50'32"E WG 12 12 
Totals 95 95 
AFLP, amplified fragment length polymorphism; WG, western group; EG, eastern group; EX, Ligia exotica group. 
*Number of individuals of which 16S rRNA is sequenced. 
>Number of individuals of which AFLP is genotyped. 
DNA is preferred to infer the overall complex phylogenetic 120E 130E 140E 
histories of the species. 
Amplified fragment length polymorphisms (AFLP) (Vos Russia 
et al., 1995) have been proven to be a useful DNA marker, 
encompassing mostly nuclear genetic variation in diverse > 42n] 
animal groups (Mueller and Wolfenbarger, 1999; Meudt and 
Clarke, 2007). Applications of these markers are effective, EGJNI i 
especially with non-model species, as there is no need for _| 
prior genomic information of the species concerned. Hence, HN 
AFLP markers are used for a variety of purposes, including 
population genetic diversity estimation, phylogeography, = 
phylogenetic analyses, genome-wide mapping, and detec- a 
tion of signatures of selection (Mueller and Wolfenbarger, 
1999; Bensch and Åkesson, 2005; Meudt and Clarke, 2007; 
Davey et al., 2011; Smith et al., 2011; Caballero et al., 
2013), and in particular, species delimitation of closely re- 140E 


lated species (Meudt and Clarke, 2007). 

To investigate whether Ligia displays similar genetic 
divergence patterns with both mitochondrial DNA and 
nuclear DNA, we evaluated the genetic variation of Ligia 
populations from South Korea, China, Japan and Russia, 
by using both mitochondrial 16S rRNA gene and AFLP 
markers. 


MATERIALS AND METHODS 


Sampling, DNA extraction, PCR and sequencing 
A total of 95 individuals of Ligia specimens were collected 
by hand from eight locations in East Asia (Table 1, Fig. 1), 
and preserved in absolute ethanol. A total genomic DNA 
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Fig. 1. Geographical locations of the Ligia sampling sites in a 
map of the East Asian region. The names of locations and relat- 
ed information are presented in Table 1. 


was extracted from the walking legs, by following Sam- 
brook and Russell (2001), or DNeasy Blood and Tissue kit 
(Qiagen, USA). Part of the mitochondrial 16S ribosomal 
RNA gene was PCR amplified, and directly sequenced using 
primers 16SAR-SF and 16SBR-SF (Palumbi et al., 1991), 
following the protocol of Jung et al. (2008). 


AFLP genotyping 
The total genomic DNA of each specimen was digested by 
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Table 2. Genetic diversity indices in sampled Ligia populations 


Genetic Variation of Sea Slaters in East Asia 


WG_EW WG_JSN WG_VB EG_UJ EG_JJ EG_JNI EG_QD EX_ZH Mean 
16S rRNA No. of substitutions 1 2 0 1 1 0 1 6 1.5 
No. of indels 0 0 0 0 0 0 0 12 1.5 
0 (7) 0.182 0.576 0 0.167 0.167 0 0.167 3 0.532 
8 (K) 0.410 0.934 0 0.393 0.393 0 0.393 2.685 0.651 
AFLP No. of polymorphic loci 380 279 201 255 214 332 235 296 449 


Proportion of polymorphic loci (%) 


AFLP, amplified fragment length polymorphism. 


using two restriction enzymes, EcoRI and Msel (Macro- 
gen, Korea) for 1 h at 37°C, followed by an additional 3 h 
at 16°C. After inactivation of the restriction enzymes, both 
EcoRI and Msel adaptors were ligated to the sticky-end of 
the restricted DNA fragments, by using T4 DNA ligase for 
1 h at 4°C. The two primers (ECO-AXX and MSE-AXX) 
were used in the pre-amplification procedure. Then, three 
types of ECO-AXX (ECO-AGG, -AAG, -AGA) with fluo- 
rescent labels (6FAM) and two types of MSE-AXX (MSE- 
ACC, -AGG) were used for the selective amplification 
procedure. The AFLP protocols used in this study followed 
Jung et al. (Jung et al., 2006, 2010; Jung, 2013) (see these 
papers for detailed information of PCR reaction mixture 
compositions and experimental conditions). The PCR am- 
plifications were performed on a GeneAmp PCR System 
2700 (Applied Biosystems, USA). The selective amplifica- 
tion products were separated in a Genetic Analyzer 3730 
(Applied Biosystems). The band sizes and genotypes were 
scored using the software GeneMarker 2.6.0 (Softgenetics, 
USA), and the existence of alleles was denoted as “1”, if the 
band intensity ranged within 100 to 60000. 


Data analyses 

Multiple alignments of 16S rRNA sequences were conduct- 
ed using Culstal W (Larkin et al., 2007), as implemented in 
the MEGA 6 packages (Tamura et al., 2013). Genetic diver- 
sity indices of each sampling sites based on 16S rRNA se- 
quences and AFLP genotypes were calculated, using ARLE- 
QUIN 3.5 (Excoffier and Lischer, 2010). The same software 
was used to infer pairwise mean Kimura 2 parameter (K2P) 
distance and Fsr among populations with permutations of 
1,000, based on 16S rRNA and AFLP, respectively. MEGA 
6 was used to depict phylogenetic trees of 16S rRNA gene 
sequences by maximum likelihood methods, with bootstrap 
values of 1,000 replicates. The best-fit models of nucleotide 
substitution for phylogenetic tree estimation were estimated 
using the same software packages. The principal coordinates 
analysis (PCoA) of AFLP genotypes was conducted in GE- 
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Fig. 2. A maximum likelihood tree among Ligia samples collect- 
ed from China, Korea and Japan based on the Tamura 3-param- 
eter obtained from 16S rRNA genes depicted with MEGA 6.06. 
The numbers near each node are bootstrap values, calculated 
using 1,000 replicates. 


NALEX 6.5 (Peakall and Smouse, 2012). 


RESULTS 


Mitochondrial 16S rRNA gene 

A total of 13 haplotypes (GenBank accession number: 
KJ802742-KJ802854) from 95 individuals belonging to the 
eight populations (Table 2). Nucleotide substitutions were 
identified at 12 sites among 497 bp, and the number of sub- 
stitutions was varied among the populations, and the indels 
with the value of 12 were only observed from L. exotica 
population (EX_ZH). Among the eight populations, WG_ 
VB and EG_JNI showed the lowest diversity values, as 0 in 
O (at) and 0 (k); whereas EX_ZH showed the highest, as 3 
and 2.685, respectively. Within the WG populations, WG_ 
JSN showed the highest diversity values. 

Three distinct genetic groups with deep divergence were 
identified in maximum likelihood trees based on the Tamura 
3-parameter, supported by bootstrap values ranging between 
95 and 100 (Fig. 2); and the Eastern group and L. exoti- 
ca were in sister relation. The pairwise tsr analysis with 
permutation of 1,000 (Table 3) between different lineages 
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Table 3. A matrix of pairwise K2P genetic distances between Ligia populations based on 16S rRNA (upper diagonal) and pairwise 
differences Fst based on AFLP data (lower diagonal), respectively 


WG_EW WG_JSN WG_VB EG_UJ EG_JJ EG_JNI EG_QD EX_ZH 
WG_EW = 0.1017 0.0083 0.9965 0.9965 0.9983 0.9965 0.9893 
WG_JSN 0.0649 z 0.1362 0.9926 0.9926 0.9942 0.9926 0.9863 
WG_VB 0.1066 0.1583 = 0.9983 0.9983 1.0000 0.9983 0.9913 
EG_UJ 0.0332 0.0366 0.1377 = 0 0 0 0.9866 
EG_JJ 0.1513 0.2070 0.1309 0.1958 = 0 0 0.9866 
EG_JNI 0.0702 0.1319 0.0648 0.0941 0.0896 = 0 0.9885 
EG_QD 0.0706 0.0059 0.1777 0.0294 0.2335 0.1484 = 0.9866 
EX_ZH 0.0113 0.0302 0.1063 0.0105 0.1762 0.0833 0.0116 = 
Significant values (p<0.05) are shaded with boldface. 
K2P, Kimura 2 parameter; AFLP, amplified fragment length polymorphism. 
ranged from 0.9916 to 1 with highly significant P, and the JEW 
values within the same lineages were 0 to 0.1364, and sta- =LE a 
tistically insignificant in all comparisons. a EX 

WG 
EG 
AFLP marker for nuclear DNA mxo a ë 
Six hundred and ninety-one loci were detected, ranging Q . "i r š = d a 
from 52 bp to 576 bp in size (Table 3), from AFLP anal- S $ +s _. * 
. vs . 
ysis of 95 individuals from eight populations. The num- ` . h z . š 
p n m k 

ber of polymorphic loci ranged from 201 (29.09%) to 380 ʻ a * sd wra Ai R 
(54.99%). Pairwise Fsr values of interspecific comparison A Po! d Slum at = ‘ 
ranging from 0.0000 (insignificant value = 0.0059) to 0.1777 . a s 6 7 a È 
were almost similar, with intraspecific values of WG (0.0649 s ‘ . 
to 0.1583) and EG (0.0000 to 0.2335) (Table 3). The PCoA r ï 
at individual level produced a two-dimensional representa- 


tion (Fig. 3), with the first two axes explaining 14.23% of 
the total variance, and among all three groups showed no 
clear cluster gathering populations. 


DISCUSSION 


The Ligia specimens found in East Asia were comprised of 
three genetically diverged groups, based on 16S rRNA gene 
differences (Table 3, Fig. 2). The mean K2P genetic distanc- 
es among the three groups (9.8-11.7) were within the range 
of average Ligia between species genetic distance (5.6- 
25.58) (Hurtado et al., 2010). The neighbor-joining tree (Fig. 
2) showed both a deep divergence between lineages, and 
highly significant bootstrap values (98-100%) at all nodes. 
These results matched well with the previous studies, using 
the same 16S rRNA gene (Jung et al., 2008; Yin et al., 2013) 
and COI gene (Yin et al., 2013) that identified the existence 
of the three distinguished genetic lineages of Ligia in East 
Asia. On the other hand, the AFLP variations between the 
populations within the groups (pairwise Fsr=0-0.23) were 
similar to the variation among the groups (pairwise Fst =0- 
0.21), in pairwise Fsr analyses. Likewise, the PCoA showed 
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Axis 1 (8.61%) 


Fig. 3. A plot representing the results of principal coordinates 
analysis (PCoA) of amplified fragment length polymorphism 
genotypes from 95 Ligia individuals, executed with GENEALEX 
6.501. Populations belonging to the Western Group, Eastern 
Group, and Ligia exotica are represented as diamond, square, 
and triangle, respectively. 


overlapping individuals belonging to the three different 
groups, with scattered patterns, with no clear clusters ob- 
served. 

This incongruence between mitochondrial and nuclear 
DNA markers may be explained by two different factors. 
First, DNAs could show different evolutionary histories, 
depending on the genomic regions of markers used. Be- 
cause of the maternal inheritance of mitochondrial DNA, 
the effective population size of mitochondrial DNA is only 
a quarter that of nuclear DNA, and displays a faster genet- 
ic divergence, compared with nuclear DNA (Avise, 2000). 
Mitochondrial DNA could depict genetic divergence oc- 
curring among the three Ligia groups; whereas there was 
not enough time to fully differentiate the nuclear DNA, and 
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these groups still share much of the ancient polymorphisms, 
in that the three groups are closely related. The distinct 
nuclear DNA genetic variation of 49 bp size difference of 
nuclear 18S rRNA gene between WG and EG (Jung et al., 
2008) supports this hypothesis, as nuclear rRNA regions 
usually experience rapid concerted evolution in metazo- 
ans (Alquezar et al., 2010). Many recent studies have also 
shown similar incongruence patterns in diverse animal 
groups, including crickets (Shaw, 2002), cichlids (Egger 
et al., 2007; Koblmiiller et al., 2007), and lizards (Leaché, 
2011). 

Second is genomic islands of speciation, which could 
explain the observed shared polymorphisms (Turner et al., 
2005). Usually, the effect of different intensity of natural 
selection pressure on loci during divergence may diversify 
the molecular evolution patterns of genes, depending on 
genomic regions. The wide spread of shared loci is due to 
only a few part of the genomic regions being highly differ- 
entiated between closely related species, these several di- 
verged genes causing a phenomenal influence, large enough 
to induce a speciation process, as speciation genes causing 
post-zygotic inviability and hybrid sterility (Coyne and Orr, 
2004; Wu and Ting, 2004; Orr, 2005; Presgraves, 2007; 
Rieseberg and Blackman, 2010; Nosil and Schluter, 2011). 

In conclusion, the analysis of 16S rRNA revealed the 
existence of the three Ligia lineages in East Asia. Despite 
the discrepancies between mitochondrial DNA and AFLP 
results, AFLP analysis showed most of the polymorphisms 
were shared among lineages. Recently, in corporation with 
the development of new generation sequencing technolo- 
gy (NGS), restriction-site associated DNA sequencing has 
been introduced as a new and alternative method for AFLP. 
Thus, the results provide a basis for future study, using NGS 
technology to scan genome-wide patterns of differentiation, 
and identify specific genes that directly influence phenotyp- 
ic changes, such as morphological, physiological, and be- 
havioural changes among lineages. 
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